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We investigate via LDA+DMFT (local density approximation combined with dynamical mean 
field theory) the manifestation of correlation effects in a wide range of binding energies in the hole- 
doped family of Fe-pnictides AFe2As2 {A = K, Rb, Cs) as well as the fictitions FrFe2As2 and a-axis 
stretched CsFe2As2. This choice of systems allows for a systematic analysis of the interplay of 
Hund’s coupling Jh and on-site Coulomb repulsion U in multi-orbital Fe-pnictides under negative 
pressure. With increasing ionic size of the alkali metal, we observe a non-trivial change in the iron 
3 d hoppings, an increase of orbitally-selective correlations and the presence of incoherent weight 
at high-binding energies that do not show the typical lower Hubbard-band behavior but rather 
characteristic features of a Hund’s metal. This is especially prominent in a-stretched CsFe2As2. We 
also find that the coherent/incoherent electronic behavior of the systems is, apart from temperature, 
strongly dependent on Jh and we provide estimates of the coherence scale T*. We discuss these 
results in the framework of reported experimental observations. 

PACS numbers: 71.15.Mb, 71.27.-|-a, 74.25.Jb, 74.70.Xa 


I. INTRODUCTION 

The nature and degree of correlations in Fe-based su¬ 
perconductors has been a subject of intensive discussions 
since the discovery of the first high-Tc iron pnictide su¬ 
perconductor in 200 ^. A significant amount of work 
has concentrated on the description of experimentally ob¬ 
served effects of correlation like larce effective masses or 
possible non-Fermi liquid behaviop^HH. Due to the multi¬ 
orbital nature of Fe-based superconductors, the Hund’s 
coupling Jh has been shown to play a key role in de¬ 
termining the behavior of these system^^Mll] However, 
there is an ongoing debate regarding the role of Jh ver¬ 
sus the on-site Coulomb repulsion U and the interpre¬ 
tation of the correla ted nature of Fe-pnictides and Fe- 
chalcogenide^^mUnHUl An important feature observed 
in many of these studies is that, depending on the elec¬ 
tronic filling, the Hund’s coupling Jh on the one hand 
renders a moderately correlated system even more corre¬ 
lated and pushes it into a bad metal regime, while on the 
other hand it can also reestablish a metallic behavi or, al- 
beit orbital selective, in a strongly correlated systenP^Hl. 
The main question for the Fe-based superconductors nar¬ 
rows down to which regime of parameters do they belong 
to and how do correlations manifest in a wide range of 
binding energies as a function of doping and/or pressure? 

Of special interest are the hole-doped AFe2As2 {A = 
K, Rb, Cs) end members of the 122 iron pnictide series. 
Compared to the parent compound BaFe2As2, the sub¬ 
stitution of Ba by K accounts for a doping of one hole 
per formula unit and it is accompanied by a complete 
supp ressio n of any structural or magnetic phase tran- 
sitiorP^li^ and by the appearance of superconductivity 


at low temperatures. This behavior is commorP^lt^ to 
all hole-doped end members AFe2As2 and they all seem 
to have a nodal gap structur ci^^ l ^^ l, which is different to 
the nodeless gap structure found in the parent BaFe2As2 
systenPMSH, Further, experimentally Bai_a;K2,Fe2As2 
is thought to undergo a coherence-incoherence transi- 
tioipadl function of temperature, that has been 
interpreted in terms of a strong increase in correla- 
tion^MEI] Measurements of the Sommerfeld coefficient 
suggest that the hole-doped end systems are among 
the stron gest correlated 122 iron-pnictide superconduc- 

t0r j^°F7 |, 

is also corroborat ed by theoretical in¬ 
vestigations on KFe2Asj '^ l ^^ l ^^ l ^^ l^. The measured Som¬ 
merfeld coefficient increases from B aFe2As 2 to KFe2As2 
by more than an order of magnitud^^ZM^^d increases 
further as K is substituted by Rb and jg view of 

these observations, the hole doped end members AFe2As2 
provide an ideal background for investigating correlation 
effects as function of negative pressure. 

In the present work we analyze via a combination 
of density functional theory (DFT) in the local density 
approximation (LDA) with dynamical mean field the¬ 
ory (LDA-I-DMFT) the electronic structure of the se¬ 
ries AFe2As2 (A = K, Rb, Cs) as well as the fictitious 
FrFe2As2 and a-axis stretched CsFe2As2 in an extended 
range of binding energies. The combination of density 
functional theory with many body methods such as dy¬ 
namical mean field theory has proven to be a very help¬ 
ful tool to investigate electronic correlation effecta^^"^. 
While a compression of the unit cell usually decreases 
correlation effects in the 122 iron pnictides^^HBl, sinii- 
larly, the opposite is to be expected when lattice param¬ 
eters are expanded since the larger interatomic distances 
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should reduce the hybridization of neighboring atomic or¬ 
bitals and lead to stronger localization of the electronic 
states. We show in this study that these considerations 
are correct only at first sight; actually, the strong corre¬ 
lation effects in these systems are mostly governed by a 
subtle interplay of Jr and U. 


II. RESULTS AND DISCUSSION 


We performed DFT(LDA) and LDA-I-DMFT calcula¬ 
tions on experimental structures of KFe2As2, RbFe2As2 
and CsFe2Asspl in tetragonal symmetry. Additionally, 
we prepared a fictitious structure for FrFe2As2 with lat¬ 
tice parameters and As height obtained by extrapolating 
the experimental parameters for AFe2As2 (A = K, Rb, 
Cs) (see Appendix A). Furthermore, we considered an a- 
axis stretched CsFe2As2 structure by slightly expanding 
the lattice of CsFe2As2 in the a-b plane to mimic an ex¬ 
perimental biaxial stretching of the crystal. This may be 
achieved experimentally by epitaxial growth of CsFe2As2 
on a substrate with slightly larger lattice parameter. The 
structural trends of this series are an expansion along a 
and c and a decrease of the As height. A detailed ex¬ 
planation of the structures and computational methods 
is given in Appendix [A} 


A. Momentum-resolved spectral function A(k, ca) 



In Fig. we show the momentum-resolved spectral 
function for KFe2As2 and CsFe2As2 as obtained within 
LDA-I-DMFT (gray density plot) at a temperature T = 
145 K and compare it to the DFT bandstructure (red). 
We estimate for both systems Fe 3 d effective masses 
to*/tolda between 2.2 and 4.1 (depending on the or¬ 
bital) which lead to a renormalization of the DFT(LDA) 
band energies and overall reduction of bandwidth. The 
average effective mass increases slightly from 2.89 to 
2 . 95 , which indicates an increase in correlation along 
the AFe2As2 series due to enhanced localization of elec¬ 
trons on Fe. Already at the Fermi level we obtain dif¬ 
fuse structures corresponding to incoherent quasiparti¬ 
cle excitations with finite lifetimes. At energies below 
— 0.5 eV all coherent features are basically washed out 
due to correlation. This effect is present in all systems 
in the series KFe2As2—>'CsFe2As2, where we observe an 
energy range of strong broadening and depletion of spec¬ 
tral weight compared to the DFT bandstructure between 
[— 2 . 0 , — 0 . 5 ] eV. At energies below — 2 eV coherent fea¬ 
tures become visible again, which correspond to the As 
p states that partially hybridize with Fe 3 d states. Even 
though the self-energy in DMFT has no momentum de¬ 
pendence, an effective momentum dependence is present 
in the results due to the momentum dependent orbital 
character of the original DFT bands. This leads to k- 
dependent broadening effects in LDA-l-DMFT. 


FIG. 1 : (Color online) Momentum-resolved spectral function 
of KFe2As2 (top) and CsFe2As2 (bottom). Strong correla¬ 
tions in these materials introduce renormalization effects as 
well as broadening in the spectral function due to finite quasi¬ 
particle lifetimes compared to DFT(LDA). We observe a no¬ 
table broadening and suppression of spectral features espe¬ 
cially in the energy range [— 2 . 0 , — 0 . 5 ] eV. 


B. Density of states at the Fermi level and 
orbital-resolved electronic filling 

We first analyze the manifestation of correlation ef¬ 
fects near the Fermi level. For that we compare in 
Fig. [| LDA-I-DMFT (Fig.[|(a)) with DFT(LDA) (Fig.[| 
(b)) orbital-resolved density of states at the Fermi level 
N[Ep) for all studied systems. The LDA-I-DMFT cal¬ 
culations show an increasing and pronounced dominance 
of 3(^22 contribution at Ep along the series. This is in 
contrast to the DFT results where the Fe 3(^22 orbital 
contribution also increases, but is much lower and only 
significant for a-axis stretched CsFe2As2 with shortest As 
height. The nonmonotonous behaviour of the Fe 
contribution in DFT is a result of the special electronic 
structure of these systems at the M point. This fea¬ 
ture is greatly suppressed by broadening effects in the 
LDA-I-DMFT calculation. 
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FIG. 2: (Color online) The density of states of the Fe 3d or¬ 
bitals at the Fermi level as obtained from (a) LDA-I-DMFT 
and as comparison with (b) DFT. The electronic correlations 
induce a marked deviation from the DFT results: The con¬ 
tribution of the Fe 3d^2 orbital is strongly enhanced, while a 
more pronounced decrease of the Fe 3dxy orbital is found to¬ 
wards the end system of stretched CsFe 2 As 2 . Also, the trend 
in the Fe orbital in DFT is completely evened out in 

LDA-bDMFT. 


FIG. 3: (Color online) Filling of the Fe 3d orbitals as obtained 
from (a) LDA-I-DMFT and (b) DFT(LDA). Correlations re¬ 
duce the overall Fe 3d filling, with the most correlated ?>dxy 
and idj.^/yz orbitals being the closest to half filling in contrast 
to the DFT result. Along the alkali metal series, the occu¬ 
pation of Fe ?)d ^2 and 3d^2_y2 orbitals reduces, while that of 
the 3dxy and 3d^^iy^ orbitals increases after a local minimum 
in RbFe 2 As 2 , which has the smallest 3dxy occupation of all 
systems studied. 


Linked to these results is the behavior of the 
LDA-I-DMFT orbitally-resolved electronic filling (see 
Fig.§. We observe that inclusion of correlation effects 
not only reduces the overall Fe 3d filling due to the addi¬ 
tional cost of the Coulomb interaction energy for doubly 
occupying a certain orbital, b ut als o introduces orbital- 
dependent effects (see section IID). While the Fe 3(1^2 
and 3dy;2_y2 filling decreases along the series, it increases 
for 3dy,zlyz and Zd^y- 


C. Spectral function A{uj) 

In order to understand the origin of the changes in 
N{Ef) and orbital filling, we show in Fig. the local 
density of states (spectral function ^(a;)) for Fe 3d^2 and 
Fe 3dj;y in the energy range [—2,2] eV for the representa¬ 
tives KFe 2 As 2 , CsFe 2 As 2 and a-axis stretched CsFe 2 As 2 . 
We chose these two orbitals since they are the most af¬ 
fected along the series: 3d22 because of the decrease in 
As height from KFe 2 As 2 to CsFe 2 As 2 and 3dxy because 
of the increase in the orbital localization with increas¬ 
ing a lattice parameter. Fe 3dz2 shows a shift of orbital 
weight to negative energies (high binding energies) and 
a narrowing of the quasiparticle-like peak structure at 
Ep from KFe 2 As 2 to CsFe 2 As 2 . Since the electronic fill¬ 


ing of the iron 3d orbitals is larger than half-filling in 
these systems, the quasiparticle-like peak is located close 
to but below the Fermi level. In Fe 3dj;2 due to the re¬ 
duction of the filling from 0.66 in KFe 2 As 2 to 0.63 in 
stretched CsFe 2 As 2 (see Fig. caused by an increase of 
electronic correlations, the quasiparticle peak-like struc¬ 
ture is shifted even closer towards the Fermi level along 
the series, which in turn leads to the observed increase of 
the density of states at the Fermi level (Fig.|^a)). Such a 
shift of the quasiparticle peak-like structure in KFe 2 As 2 
was already noted in Ref. [Ml For the 3dxy orbital, the 
quasiparticle peak is much closer to the Fermi level since 
the filling is closer to half filling at n^y « 0.59. This 
orbital shows a strong suppression of the quasiparticle 
peak (up to « 30%) from KFe 2 As 2 to the a-axis stretched 
CsFe 2 As 2 which points to an important increase of de¬ 
coherence along the series. This increase of decoherence 
will be studied in more detail in the next section. Reduc¬ 
tion of the maximum of the quasiparticle peak combined 
with a slight change of its position results in the almost 
constant density of states at the Fermi level for the 3dxy 
orbital observed in Fig. |(a). 

Additionally, a shoulder-like feature appears in the 
spectral function at 1 eV and — 1 eV in CsFe 2 As 2 , closely 
resembling the typical spectral function shape of a quasi¬ 
particle peak and a lower and upper Hubbard band. 
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FIG. 4: (Color online) The density of states for (a) the Fe 
3d^2 orbital and (b) the Fe Sd^y orbital as obtained from 
LDA-I-DMFT for the three compounds KFe 2 As 2 , CsFe 2 As 2 
and an a-axis stretched CsFe 2 As 2 . Note the emergence of a 
Hubbard-like peaks around —1.2 eV and -|-1 eV for 3d22 and 
around —1 eV and -1-1 eV for Sd^y. With arrows we mark the 
trend of the changes in the spectral function along the series. 


These features do not correspond to any property found 
in the non-interacting DOS and are purely an effect of 
correlations and, at first sight, are similar to the emer¬ 
gence of Hubbard bands as a function of U in strongly 
correlated systems. This is also in agreement with the 
3dxy orbital being the strongest correlated one, whereas 
these features are far less developed in the 3^22 and other 
orbitals (see Appendix [D|). However, as we will show in 
section m these peaks do not behave as expected for 
Hubbard bands in a one-band Hubbard model. 


D. Effective masses and quasiparticle lifetimes 

In order to quantify the change in correlation along the 
series 24Fe2As2 we plot in Fig. the orbitally resolved 
quasiparticle lifetimes and effective masses . The 
effective masses for Fe 3^22 and 3dx2-y2 increase along 
the 24Fe2As2 series but remain constant or even slightly 
decrease for ^d^zjyz and Sd^y This last result cannot 
be explained solely by the behavior of the DFT-derived 
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FIG. 5: (Golor online) (a) Quasiparticle lifetimes given by 
Tm = —{ZmlrDSrn{iO'^))~^ of the Fe 3d orbitals along the 
AFe 2 As 2 series and (b) the corresponding effective masses 

m* 

^LDA ■ 


tight binding parameters (see Appendix and . 

Information on the origin of this behavior can be ob¬ 
tained from the quasiparticle lifetimes which strongly de¬ 
crease for all orbitals along the series (Fig.[^(a)). This in¬ 
dicates that the coherent quasiparticle picture, being the 
basis for the calculation of the effective masses, becomes 
less appropriate along the AFe 2 As 2 series, and accord¬ 
ingly, the effective masses obtained by this procedure are 
an underestimation of the true value. This result shows 
that already at the temperature of T = 145 K these sys¬ 
tems are quite incoherent. With increasing lattice param¬ 
eter incoherence significantly increases, albeit leading to 
an orbital dependent change in localization of the Fe 3d 
electrons. In particular, we find a pronounced decrease 
of the quasiparticle lifetime from KFe 2 As 2 to RbFe 2 As 2 , 
which we attribute to the competing effects of decrease of 
Fe-Fe direct hopping and increase in Fe-Fe indirect hop¬ 
ping through As, which seem to have a crossover point 
between RbFe 2 As 2 and CsFe 2 As 2 (see Appendix [B|) . 

In Fig. we show the temperature dependence of 
the quasiparticle scattering rates and effective masses of 
KFe 2 As 2 . We find a similar temperature dependence as 
also observed for multi-orbital SrRuOa and CaRuOa in 
Ref. |46l As the temperature is lowered in the calcu¬ 
lation, the quasiparticle picture, which is suppressed at 
high temperatures, is partially restored, leading to an 
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FIG. 6: (Color online) Temperature dependence of the (a) 
scattering rates of the Fe 3d orbitals of 

KFe 2 As 2 . The shaded area indicates the coherent domain. 
The coherence temperature estimate T* as deduced from the 
scattering rate is quite low, located around 50 K. (b) The 
corresponding effective masses. 


exponential increase in the quasiparticle lifetimes. Still, 
even at the lowest studied temperature of T = 72 K 
the width (scattering rate) of the quasiparticle peak 
— ZmImEm(*0“*') is still larger than the temperature, in¬ 
dicating that coherent quasiparticle excitations are still 
in the minority. From our results for the scattering rate 
we estimate that the temperature T* for the incoherent- 
coherent phase transition for our chosen interaction pa¬ 
rameters C/ = 4 eV, Jh = 0.8 eV is around 50 K, and is 
even lower for RbFe 2 As 2 and CsFe 2 As 2 . This is in qual¬ 
itative agreement with previous LDA-I-DMFT estimates 
for KFe 2 AsJ^and with magnetic susceptibility and ther¬ 
mal expansion measurements which pr edict th e coher¬ 
ence scale of KFe 2 As 2 to be around lOOKP^IsniSl and even 
lower for RbFe 2 As 2 and CsFe 2 As J^°*^^ l We expect that 
inclusion of the full rotationally invariant Hund’s cou¬ 
pling, which is beyond the scope of the present work, will 
shi ft the t ransition to higher temperatures in the calcula- 
tioipEBll and, therefore, closer to experiment. The con¬ 
sequences of inclusion of rotationally invariant interac¬ 
tions have been extensively discussed in previous works in 
the framework of the two-bancP^ and three-bancP^ Hub¬ 
bard models. The computational effort of including these 
terms in the five-band cases studied here is beyond the 
scope of the present work. However, we performed some 
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FIG. 7: Sketch of some of the most likely Fe atomic orbital 
configurations in hole-doped 122 iron pnictide systems. The 
electronic filling is indicated by n and the total spin by S. 


benchmarking calculations with fully rotational Hund’s 
coupling and find that the coherence transitio n shift s to 
higher temperatures as found in other studied^*^*^. In 
this sense the presented results can be considered as a 
lower bound to the experimental observations. 

The increase of effective masses at lower tempera¬ 
tures in our calculation is precisely the effect of restoring 
the coherent quasiparticle picture, so that the effective 
masses at the lowest temperature (72 K) we investigated 
can be considered as the closest approximation to the 
true values, i.e. 6.1 {3dxy), 4.6 (id^z/yz), 3.4 (3c?^2) and 
2.5 {3dx2—y2). 

Combining these observations along the series we con¬ 
clude that alkali 122 systems show typical signs of strong 
correlations but, in the studied range of temperatures, 
are actually quite deep in the incoherent bad metal region 
with a well defined, albeit strongly suppressed, quasipar¬ 
ticle peak. 


E. Nature of the Fe 3d wavefunction 


We proceed now with an analysis of the wave function 
in terms of the Fe atomic basis states, similar to what 
has been done for other pnictide systemd^. In Fig. we 
show a sketch for a few typical orbital configurations in 
these systems with the orbital splittings obtained from 
the downfolded charge self-consistent LDA Hamiltonian. 
Since the nominal electronic filling for Fe 3d for these 
systems is 5.5 electrons per Fe, one would expect atomic 
Fe 3d states with 5 and 6 electrons to be the most likely 
states. However, as shown in Fig.[^ the actual average Fe 
3d Hlling is slightly larger than 6 and this can be analyzed 
in the histogram of Fig. . For the interacting system 
there is a non-trivial competition between the energy con¬ 
tribution due to the crystal field splitting, which prefers 
to occupy the lowest states first, the on-site U interaction, 
which tends to decrease the filling of the localized states 
and the Hund’s coupling Jh, which prefers orbital states 
with maximum total spin. When the Hund’s coupling Jh 
is large compared to the total crystal field splitting, the 
high-spin states will have the highest probability and the 
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FIG. 8: (Color online) The histogram of the Fe 3d atomic 
state for KFe 2 As 2 at T = 145 K. The probability corresponds 
to the fraction of time the Fe 3d orbitals spend in one of the 
2^® = 1024 possible states. Within the interval of constant 
electron number N the states are sorted by increasing energy, 
i.e. the leftmost states within an interval correspond to high- 
spin states, while the rightmost states correspond to low-spin 
states. 
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low-spin states will be suppressed. This is indeed true 
for the hole-doped 122 iron-pnictides. 

In Fig. 1^ we show for KFe 2 As 2 the atomic histogram 
of the Fe 3d shell, i.e. the projection of the wave func¬ 
tion onto the Fe 3d atomic basis states. The states are 
sorted by the number of electrons N, and inside the 
interval of constant filling the states are sorted by en¬ 
ergy. Because of the Hund’s coupling, the leftmost states 
in such an interval correspond to the high-spin states, 
while the rightmost states correspond to low-spin states. 
The probability assigned to each state corresponds to the 
fraction of time in the calculation the 3d orbitals spend 
in a specific configuration. Due to the Hund’s coupling 
the high-spin states clearly dominate the histogram, al¬ 
though their probabilities are quite low with less than 
2%. Even for the low-spin states at higher energy the 
probability, while being smaller, is never close to zero, 
which causes the Fe 3d orbitals to visit a large number 
of accessible states over time even when they are much 
higher in energy. Since all these states contribute with a 
finite fraction to the Fe 3d spectral function, no well de¬ 
fined atomic like excitations can be expected. This leads 
to the observed suppression of the quasiparticle peak and 
subtle high-binding energy features in the spectral func¬ 
tion without well defined Hubbard band s, which is the 
distinctive property of Hund’s metal^^^^. 

This behavior becomes even more pronounced when 
we perform this analysis along the AFe 2 As 2 series. In 
Fig.[^ (a) we analyze the most likely atomic states from 
the histogram of the Fe 3d orbitals for AFe 2 As 2 . For 


FIG. 9: (Golor online) (a) Probabilities of the most likely 
atomic states of the Fe 3d atomic orbitals for the AFe 2 As 2 
series. The probability corresponds to the fraction of time the 
atom in the DMFT calculation spends in a specific state, (b) 
The summed probability of all atomic states with a given total 
spin S. High-spin states become more likely for increasing 
lattice parameter a, while the probability of low-spin states is 
reduced. 


all systems, six out of the seven most likely atomic states 
are solely composed of the maximal high-spin states with 
S = 2.5 and S = 2. For the earth alkali (undoped) 122 
iron-pnictides like BaFe 2 As 2 the atomic ground state of 
the Fe atom with a valence charge of 6 has a maximum 
possible spin of S' = 2. Since for the systems studied 
here one electron per formula unit has been removed by 
hole doping, the probability for the fully polarized half- 
filled S = 2.5 state with 5 electrons in the Fe 3d or¬ 
bitals is among the most likely states with a comparably 
high probability of 1.4% (or 2.8% when accounting for 
spin degeneracy) in KFe 2 As 2 , that increases up to 1.8% 
(or 3.6%) in the stretched CsFe 2 As 2 system. Along the 
series the low-spin states become suppressed while the 
high-spin states increase in probability, as can be seen in 
Fig.|^(b). 

Generally, the systems can reduce their energy by as¬ 
signing a higher probability to high-spin states due to 
the Hund’s coupling Jh- This leads to a significant 
increase of localization caused by the orbital blocking 
mechanisirP. Since Jh enforces a high-spin state, or- 
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bital mixing is greatly suppressed compared to a van¬ 
ishing Hund’s coupling where high- and low-spin states 
would have equal energies and, therefore, probabilities. 
This is the typical behavior of a so-called Hund’s metal, 
in which the electronic correlations are much more sen¬ 
sitive to the value of Jh than to the on-site Coulomb 
term U. Therefore, in the hole-doped end systems like 
CsFe 2 As 2 and especially the o-axis stretched CsFe 2 As 2 
the Hund’s coupling becomes the most important inter¬ 
action that governs the physical properties of these sys¬ 
tems. 


density of states. 

Our results confirm the general picture of the iron pnic- 
tides being ’’Hund’s metals” with strong orbital separa¬ 
tion, especially for the strongly correlated hole-doped end 
systems considered in this study. In this case, a slight in¬ 
crease of Jh renders the system much more incoherent 
and ’’bad metal”-like for the same value of U, while the 
spectral weight at the Fermi level differs strongly between 
the orbitals but remains finite even for larger values of 
U. 


F. Dependence on U and Jh 


III. CONCLUSIONS 


In order to investigate the effects of U and Jh more ex¬ 
plicitly and to determine the nature of the peak/shoulder 
at [—1.5,—!] eV we performed calculations for different 
interaction parameters for the most correlated case, the 
a-axis stretched CsFe 2 As 2 system. For computational 
efficiency these calculations were done at higher temper¬ 
ature j3 = AQ eV“^. While the height of the quasiparticle 
peak is reduced at higher temperatures, the behavior of 
the spectral function at [-2,-1] eV is quite robust. We 
considered on-site Coulomb values C/ = 4, 5 and 6 eV and 
Hund’s couplings of Jh = 0.8, 1.0 and 1.2 eV. In Fig. [T0| 
and 11 we show the spectral function A{u!) for Fe 8^22 
and 3dxy An increase of U from 4 eV to 6 eV at a fixed 
Jh implies only moderate changes in the spectral func¬ 
tion in general for all Fe 3d orbitals. The Hubbard-like 
shoulder at —1.2 eV becomes more pronounced for larger 
U values and its maximum moves only very slightly to 
negative energies (—1.4 eV). Due to particle-hole asym¬ 
metry, we obtain a quasiparticle like peak slightly shifted 
away from the Fermi level. On the other hand, an in¬ 
crease in the Hund’s coupling Jh for fixed U immedi¬ 
ately renders the system very incoherent, with a strong 
increase in the scattering rate and a reduction of the 
quasiparticle lifetime, leading to a strong suppression of 
DOS at the Fermi level and a significant shift of spectral 
weight to lower energies, forming a broad lower Hubbard 
band located between —4 and —5 eV. Finally, the com¬ 
bined effect of U and Jh yields an even more well defined 
lower Hubbard-like band at around —4.5 eV. While this 
characteristic dependence on U and Jh is very similar 
for all Fe 3d orbitals, we still find a strong orbital selec¬ 
tion regarding the remaining spectral weight at the Fermi 
level. Especially the Fe ‘Jtdxy orbital is almost gapped at 
the Fermi level for the largest interaction values consid¬ 
ered (see Fig. 11), while the least correlated Fe Zdx'^-y^ 
orbital experiences basically no suppression of the DOS 
at the Fermi level regardless of the interaction parame¬ 
ters and instead retains a well defined quasiparticle peak 
for higher values of U and Jh (see Appendix]^. We also 
checked the case of negligible Hund’s coupling by setting 
Jh = 0 , which recovered the coherence properties even 
at r = 300 K, with low effective masses around 1.4 and 
a spectral function that resembled quite well the DFT 


From our analysis of the electronic properties within 
LDA-I-DMFT in a wide range of binding energies, we con¬ 
clude that along the isoelectronic doping series AFe 2 As 2 
{A = K, Rb, Cs) as well as the fictitious FrFe 2 As 2 and 
a-axis stretched CsFe 2 As 2 , which shows a monotonous in¬ 
crease of the a lattice parameters and a decrease of the As 
z height, correlation and incoherence of the Fe 3d orbitals 
increase, albeit orbitally selective, and the systems show 
clear features of a Hund’s metal. In this case the Hund’s 
coupling plays the major role and renders these materials 
much more incoherent than expected from the value of 
the Coulomb repulsion U alone. While the most corre¬ 
lated orbitals {dxy) show features that resemble those of 
being close to an orbital selective Mott transition, spe¬ 
cially for a-stretched CsFe 2 As 2 , the system is quite deep 
in the incoherent bad metal regime with a finite spec¬ 
tral weight at the Fermi level even for C7 = 6 eV and 
Jh = 1.2 eV. Experimentally, we predict that an increase 
of the Fe-Fe distance in CsFe 2 As 2 by stretching will in¬ 
duce an orbital dependent increase in correlations and 
incoherence of the Fe 3d orbitals, where the Fe 3dz'i and 
Fe 3dxy orbitals are strongly but not fully localized and 
the other Fe 3d orbitals retain a bad metallic behaviour. 
From our results we estimate the coherence temperature 
to be located around 50 K in KFe 2 As 2 and even lower for 
RbFe 2 As 2 and CsFe 2 As 2 in qualitative agreement with 
the experimental observations. These features make the 
hole doped end systems of the 122 iron pnictides, namely 
KFe 2 As 2 , RbFe 2 As 2 and especially CsFe 2 As 2 and a-axis 
stretched CsFe 2 As 2 a valuable test bed to study the be¬ 
havior of strongly correlated Hund’s metals and orbital- 
selective bad metallicity and its interplay with supercon¬ 
ductivity. 
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FIG. 10: (Color online) The density of states for the Fe 3d^2 orbital of the stretched CsFe 2 As 2 compound as a function of 
the on-site Coulomb repulsion U and Hund’s coupling Jh- (a) An increase only in U leads to an increase in renormalization, 
i.e. effective masses and a pronounced Hubbard-like peak at —1.5 eV but no other qualitative changes are observed, (b) The 
Hund’s coupling Jh greatly increases the decoherence of the electronic states at the Fermi level and leads to a significant shift 
of spectral weight down to lower energies, (c) The combined effect of U and Jh is qualitatively very similar to an increase in 
Jh alone. 
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FIG. 11: (Color online) The density of states for the Fe orbital of the stretched CsFe 2 As 2 compound as a function of the 
on-site Coulomb repulsion U and Hund’s coupling Jh- Similar to the Fe 3d^2 orbital, we see that the shift of spectral weight 
to lower energies is almost exclusively dependent on the Hund’s coupling Jh- Being the most correlated orbital, the DOS of 
the Fe Sd^y orbital at the Fermi level is almost gapped for 1/ = 6 eV and Jh = 1.2 eV. 


Appendix A: METHODS 

For our fully charge self-consistent LDA-FDMFT cal¬ 
culations (see Ref. |33] for detailed explanation) we use 
the structural parameters from Ref. 0^ for the tetragonal 
structures of KFe 2 As 2 , RbFe 2 As 2 and CsFe 2 As 2 at room 
temperature. Due to the almost perfect linear depen¬ 
dence of the lattice parameters as a function of atomic ra¬ 
dius, we further use linear extrapolation to obtain struc¬ 
tural parameters for fictitious FrFe 2 As 2 , avoiding pos¬ 
sible ambiguities from DFT-based relaxation methods 
which do not work satisfactorily for these systems. Ad¬ 
ditionally, we prepare a structure for CsFe 2 As 2 that is 
extended along the a/&-axis by 3% and has a reduced 
relative As z height of 2% to mimic a small expansion of 
the lattice. The expansion is performed in both the x- 
and y-direction so that tetragonal symmetry is preserved. 


The lattice parameters and As z position are shown in 

Fig.[g 

For the DFT calculations we used the WIEN2 k!^ im¬ 
plementation of the full-potential linear augmented plane 
wave (FLAPW) method in the local density approxima¬ 
tion. The Kohn-Sham equations were solved on 726 k- 
points in the irreducible Brillouin zone, resulting in a 
21 X 21 X 21 fc mesh in the conventional Brillouin zone. 
A local orbital basis was obtained by a projection of 
the Bloch wave functions to the localized Fe 3d orbitals, 
using our implementation of the method described in 
Refs. l9l5Tl Please note that a coordinate system which is 
rotated by 45° around the z-axis with respect to the con¬ 
ventional lA/mmm unit cell is used. Thus, the x- and 
y-axis point towards neighboring Fe atoms. The energy 
window for the bands to be considered for projection was 
chosen to be [—6,13] eV, with the lower boundary lying 
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FIG. 12: (Color online) The structural parameters for the 
considered materials KFe 2 As 2 (K), RbFe 2 As 2 (Rb), CsFe 2 As 2 
(Cs) (from Ref. 1451) . and the fictitious systems FrFe 2 As 2 (Fr) 
and stretched CsFe 2 As 2 {Cs stretch) ■ (a) The a-parameter and 
absolute height of the As-atom above the Fe plane, (b) The 
c-lattice parameter. 
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FIG. 13: (Color online) Eight largest effective Fe-Fe hop¬ 
ping parameters t obtained from a 10-band tight binding fit. 
Also shown are some second nearest neighbor parameters in¬ 
dicated by the label NN2. The non-trivial evolution of the 
hopping values with increasing lattice parameter and reduced 
Asz height leads to different degrees of localization and effects 
of correlation in the Fe 3d orbitals along the series KFe 2 As 2 , 
RbFe 2 As 2 , CsFe 2 As 2 and FrFe 2 As 2 . 


in a gap in the density of states (DOS). Consequently, 
35 bands on average were taken into account for the pro¬ 
jection, resulting in a representation of the fc-dependent 
and local non-interacting spectral function for each or¬ 
bital that is indistinguishable from the DFT result in the 
chosen energy window. 

The DMFT impurity problem was solved with the 
continuous-time quantum Monte Carlo method in th e hy- 
bridization expansiorP^ as implemented in the 
project. In the calculations we used an inverse temper¬ 
ature of ,5 = 80eV“^, corresponding to the temperature 
of 145 K, unless stated differently. A total number of 
at least 50 x 10® Monte-Carlo sweeps were performed for 
each solution of the impurity model and up to 90 x 10® 
sweeps for the larger interaction parameters. For the 
double count ing correction we used the nominal double 
countin j®® ! ®® !, which has been shown to yield significantly 
better agreement with photoemission experiments^, es¬ 
pecially for low and high binding energies, while other 
methods like the FLlP^EH double counting scheme over¬ 
estimate the valence charge and underestimate a possible 
Mott gapP^. The interaction parameters were used in the 
definition of the Slater integralP^ with U = and 
Jh = {F^ + F^)/14. For the on-site Coulomb interac¬ 
tion we considered a value of 17 = 4eV and for Hund’s 
rule coupling Jh = 0.8 eV, unless stated differently. We 
calculate the effective masses directly from the impurity 
self-energy via 


m 


= 1 - 


9ImE(fa;) 


wlda 


duj 


(Al) 


>0+ 


with the quasiparticle weight being defined as the inverse 


of the effective mass Z = 


m* 


-1 


The continuation 


of the Monte Carlo data to the real axis was done by 
stochastic analytic continuatiorP^I. 


For obtaining the hopping matrix elements we ob¬ 
tained a tight-binding Hamiltonian from projective Wan- 
nier functionPi^ from DFT, generated by the all-electron 
full-potential local orbital (FPLO))^ code, using a 10 
(16) orbital model, including the Fe 3d only (10-orbital 
model)P^or Fe 3d and As 4p orbitals (16-orbital model). 


Appendix B: 10-band tight binding model 

In order to quantify at the level of DFT the effects 
of negative pressure introduced by isovalent doping in 
AFe 2 As 2 {A = K, Rb, Cs, Fr), we calculated the Fe-Fe 
hopping matrix elements via projective Wannier func¬ 
tions. The absolute values of the hopping parameters 
are plotted in Fig. |13[ There are two main contributions 
that affect the values of the Fe-Fe hopping parameters: 
First, the increase of the interatomic distances due to in¬ 
creasing atomic radii of the alkali ions implies a decrease 
of the direct Fe-Fe hopping. Second, due to elongation of 
the Fe-As tetrahedron, the As atom moves closer to the 
Fe-Fe plane. This reduction of the As^ height leads to 
an increase of the indirect hopping along the path Fe-As- 
Fe. The total contribution of these two effects translates 
into a non-trivial behavior of the Fe-Fe effective hoppings 
along the doping series (K,Rb,Cs)Fe 2 As 2 . The Fe 3dxy- 
Fe 3dxy effective hopping is the smallest in KFe 2 As 2 due 
to the almost perfect cancellation of the two contribu¬ 
tions as pointed out in Ref. 01 As the lattice parame- 
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fe^-3cl^_yf ->Fe^- 3d^_y^ 
Fe'-3d/_yS ->Fe^- 3d/ 
Fe^-3d^^^^->Fe2- 3d^y^ 
Fe'-3dxy->Fe^- 3d^y 
Fe^-3d^y->Fe^- 3d^y^ 
Fe'-3d/->Fe^- 3d/ 


FIG. 14: (Color online) The first six largest hopping param¬ 
eters t for the Fe-Fe hopping obtained from a 16-band tight 
binding fit. Fe“ denotes the i-th atom out of the two equiva¬ 
lent iron atoms in the irreducible Brillouin zone. 


ter increases, the indirect Fe-As-Fe hopping contribution 
outweighs the contribution coming from the direct Fe- 
Fe orbital overlap. This causes a slight increase of the 
hopping parameters from KFe 2 As 2 to CsFe 2 As 2 . These 
two contributions are very similar for hoppings between 
Fe id^z/yz orbitals. The trend in the hoppings between 
Fe and 2>dx2_y2 orbitals is less affected by changes 
of the indirect hopping contribution and shows a small 
overall decrease in the hopping to the neighboring Fe Zd 
orbitals. 


Appendix C: 16-band tight binding model 


and a-axis stretched CsFe 2 As 2 . Compared to the other 
Fe 3d orbitals they are less affected by an increase of 
the lattice parameter a. Similar to the 3d^2 and id^y 
orbital a small Hubbard-like peak becomes more pro¬ 
nounced in the id^z/yz orbital, while the 3dx2-y2 orbital 



-2 -1.5 -1 -0.5 0 0.5 1 1.5 2 

( 0 [eV] 


In Fig. we show the six largest Fe-Fe hopping pa¬ 
rameters obtained from a 16-band tight binding fit, en¬ 
compassing the Fe 3d and As 4p orbitals. The overall 
monotonous decrease resembles the increase of the in¬ 
teratomic distance that leads to a reduced overlap of 
the neighboring Fe 3d orbitals. As noted in the main 
text, the indirect hopping through the As 4p orbitals has 
an important effect on the effective hopping parameters. 
Taking only the direct Fe-Fe hopping into account, we 
observe the expected decrease of the hopping parameters 
which resembles the reduced hybridization as the lattice 
parameters are increased. In combination with the indi¬ 
rect hopping via the As atom, this leads to a nontrivial 
behaviour of the effective hopping parameters. 


Appendix D: DOS for other orbitals 


In Fig. 15 we show the local spectral function of the 
Fe 3dxz/yz and idx 2 _y 2 orbital for KFe 2 As 2 , CsFe 2 As 2 


FIG. 15: (Golor online) The density of states for (a) the Fe 
?>dxy orbital and (b) the Fe 3d^2_y2 orbital for [7 = 4 eV and 
Jh = 0.8 eV as obtained from LDA-fDMFT for the three 
compounds KFe 2 As 2 , GsFe 2 As 2 and an expanded structure 
of GsFe 2 As 2 . 


shows the opposite trend, increasing its spectral function 
at the Fermi level at the cost of decreasing it at negative 
energies. 

In Fig. 16 we show the dependence of the Fe Sd^ziyz 
and 3dx2-y2 orbital spectral function on U and Jh- The 
results are very similar to the other orbitals, with the ef¬ 
fect of increasing U being much less extreme than that of 
Jh- While at higher U the spectral functions still mimics 
the LDA result, an increase in Jh results in a large shift 
of spectral weight to negative energies. The 3d^2_y2 or¬ 
bital shows the smallest degree of correlations and is the 
only orbital that retains a well defined quasiparticle peak 
at the Fermi level even for larger interaction parameters. 
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FIG. 16: (Color online) The density of states for the Fe and 3d^2_y2 orbital of the stretched CsFe 2 As 2 compound as 

a function of the on-site Coulomb repulsion U and Hund’s coupling Jh- (a)-|-(d) An increase only in U leads to a slightly 
better pronounced lower Hubbard band-like feature but otherwise no qualitative changes to the high energy features, (b)-l-(e) 
The Hund’s coupling Jh greatly increases the decoherence of the electronic states at the Fermi level and leads to a signihcant 
shift of spectral weight down to negative energies, (c)-l-(f) The combined effect of U and Jh is qualitatively very similar to an 
increase in Jh alone. 
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